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Abstract 

The rebound indentation test consisting of the displacement-controlled and load- 
controlled stages is considered for a frictionless cylindrical flat-ended indenter. The 
mechanical behavior of an articular cartilage layer sample is modeled in the frame- 
work of viscoelastic model with time-independent Poisson's ratio. Closed-form an- 
alytical expressions for the contact force (in the displacement controlled stage) and 
for the indentation displacement (in the load-controlled stage) are presented for an 
arbitrary viscoelastic solid model. The case of standard viscoelastic solid model is 
considered in detail. It has been established that the rebound displacement (that is 
the indentation displacement in the load-controlled stage) does not depend on the 
relaxed elastic modulus and Poisson's ratio as well as on the layer's thickness. 
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1 Introduction 



Indentation tests have been used for decades to identify mechanical properties of articular 



cartilage (Hayes et al. , 1972) and to assess its viability (Broom and Flachsmann, 2003) 



There is a wide variety of indentations tests and their possible applications. First of all, 
indentation tests can be classified with respect to the type of indenter used (cylindrical 
fiat-ended, spherical, etc.). In the simplest case of cylindrical indenter, the contact area 
is supposed to remain constant during the indentation process, and only two mechanical 
variables, namely, the contact force, -P(t), and the indenter displacement, w{t), can be 
recorded as functions of time during the depthsensing indentation. Second, indentation 
tests can be viewed according to the main mechanical parameter determined (e. g., short- 



time bulk modulus (Hori and Mockros, 1976), indentation stiffness (Lyyra et al. , 1999) 



pulsatile dynamic stiffness (Scandiucci de Freitas et al. , 2006)) that, of course, depends 
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on the mathematical model employed. For instance, the derivative dP/dw is known as 
the indentation stiffness, and will be a variable quantity for time-dependent materials like 
articular cartilage. Third, indentation tests can be classified according to protocol types 
(e.g., multiple ramp-and-hold protocol for creep test (Oyen, 2005)). This means that the 



time-history of loading is important. Finally, depending on the primary variable parameter 
chosen, one can distinguish displacement-controlled or load-controlled indentation tests, 
when wit) or P{t) is assumed to be controlled, respectively. In particular, linear and non- 
linear monotonic displacement-controlled and force-controlled indentation loading paths 



were considered recently by Cheng and Yang (2009). 



In their recent paper. Brown et al. (2009) employed a hybrid-type test consisting of two 



stages, one of which is displacement-controlled, while the other can be regarded as load- 
controlled. In the present note, we apply the viscoelastic layer model for evaluating the 
so-called rebound indentation test and clarifying its applicability for assessing viability of 
articular cartilage. 



2 Rebound indentation test 



We consider a cylindrical, fiat-ended indentation test, which is composed of two stages. In 
the first stage, called the indentation phase, the sample is subject to loading at a constant 
speed vq to achieve an indentation depth of wq. That is the indenter displacement is 
assumed to be specified according to the law 



w 



(1) 



(t) = VqI, 0<t<tr, 



(1) 



The maximum indenter displacement at the end of the first stage is given by 

Wq = Votm- (2) 



Hence, specifying the values of vq and wq, we obtain the duration of the first stage tr, 
wq/vq. Note that instead of an imposed maximum displacement. 



Wo, 



Brown et al. (2009) 



operated with a prescribed value of maximum strain (deformation) Wo/h, which depends 
on the layer sample thickness h. 

Further, we assume that at the indentation depth wo the load is immediately removed 
and the second stage, called the recovery phase, lasts for a theoretically indefinite time. 
In the recovery phase, we have 



(3) 



We distinguish in the notation the displacements, '{t) and '{t), and the contact 
loads, P^^\t) and P^'^\t), corresponding to aforementioned two test stages. In the first 
stage, the function w^^\t) is specified by Eq. ([l|, while P'-^^t) is unknown. On the con- 
trary, the displacement function w^'^\t) is unknown, whereas the contact load P^'^\t) is 
specified by Eq. ([s]), in the second stage. 



3 Evaluation of the rebound test using a viscoelastic model 



Let us assume that an articular cartilage sample may be modeled as a viscoelastic layer of 
thickness h bonded to a rigid substrate. For the sake of simplicity, we neglect friction and 
assume that Poisson's ratio, z/, of the layer material is time independent. These are stan- 
dard assumptions made in most of the mathematical models. Then, based on the solution 



to the elastic contact problem of frictionless indentation of an elastic layer (Lebedev and 



Ufliand, 1958 Hayes et al. , 1972) and applying the elastic- viscoelastic correspondence 



principle (Lee and Radok, 1960), one can arrive at the following equation between the 



indenter displacement w{t) and the applied force P{t): 



2a r 
P{t) = ^— ^/t(a,z/) J E{t 



r)^(r)dr. 



(4) 



Here, a is the contact radius, v) is a dimensionless factor accounting the thickness 
effect through the ratio a = a/h ( [Hayes et al. 1972), E{t) is the relaxation modulus, t is 
the time variable, t = 0— is the time moment just preceding the initial moment of contact. 
It has been tacitly assumed that w{t) = for t < 0. Note that Eq. ^ was previously 



used in a number of studies (Zhang et al. , 2004 Cao et al. , 2009). 



Let Eoo be the relaxed elastic modulus, that is the limit of modulus E{t) at t — oo. 
Taking into account the relation E{t) = Eoo'^lt), where \E'(t) is the relaxation function, 
we rewrite Eq. (|4]) in the following form: 



Pit) 



2aE^ 



r)^(r)cir. 



(5) 



Considering Eq. ([s]) as an integral equation with respect to w{t), one arrives at the inverse 
relationship 



l-u^ 1 I- .dP 



w(t) = — . r / $(t-r) — (r)rfr, (6) 



0- 



where is the creep function. With the appropriate relaxation function \l'(t) and creep 
function Eqs. ^ and (|6| can be used in the general case of viscoelastic solid. 

In the indentation phase, according to Eqs. ([T]) and ([s]), we will have 

p{l)^t)='^^^^CC,v)vo j ^{t-T)dT, 0<t<tm. (7) 





The maximum contact load Pm = P^^\tm) is given by 



Pm = v)v, J ^{t) dr. (8) 



In the recovery phase, in accordance with Eq. (|6]), we can write 



In view of ([s]), the last integral in ^ is equal to —Pm^{t — tm), where Pm is defined 
by ([s]). The derivative dP^^^dr entering the first integral in ^ can be computed by 
differentiating the both sides of Eq. ([T]). Thus, taking also into account Eq. ([s]), we can 
rewrite Eq. ^ as follows: 

w^^\t) =Voj m - r) - - t„)]^(r) dr. (10) 



It is interesting and very important to observe that, in view of (10), the rebound displace- 
ment w^'^\t) does not depend on the layer material constants E^o and v as well as on the 
layer's thickness h. 



4 Example. Standard viscoelastic solid model 



If the layer's material follows a linear standard three-parameter viscoelastic solid model, 
we will have 



Fit) = E^{1 - (1 - 1/p) exp(-t/(pr,))}. 



11) 



where Tg is the characteristic relaxation time of strain under applied step of stress, p is 
the ratio of the relaxed elastic modulus Eoo (the limit of modulus E(t) at t — > oo) to the 
unrelaxed elastic modulus Eq (modulus E{t) at t = 0), i. e., p = Eoo/Eq < 1. 



According to Eq. (11), we have 



vl/(t) = 1 - (1 - 1/p) exp(-t/(pr,)), $(t) = 1 - (1 - p) exp(-t/r,). 



(12) 



After the substitution of the expressions ( pL2| ), Eq. ( [10| yields 

^'-]<tm- pTs + pTs exp( - 



w^'^\t) = t>o(l — p) exp(^- 



pTs + pr^expf-— ) 1. 



(13) 



It is interesting and important to point out that in view of ( 13 ) the rebound displacement 
w^'^\t) depends on the time variable t only through the exponent factor exp(— (t — tm)/''"s)- 



5 Rebound strain. Comparison with experimental observations 



Following Brown et al. (2009), we determine the rebound strain, 5r(^); as the relative 
elastic reverse displacement that the viscoelastic layer recoves to from the moment of 
unloading. In other words, we put 



Srit) = ^, t > tr, 



(14) 



Now, taking into account Eqs. ([!]) and ([2]), we rewrite Eq. (14) in the form 

«;(2)(t) 



Brit) = no 



h 



(15) 



where uq = wo/h is the maximum strain at the end of the loading stage, and it is assumed 
that t >tjn. 
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Figure 1 . Rebound versus time behavior before and after one hour trypsin treatment following 
0.1 s~^ indentation. Experimental data were obtained by Brown et al. (2009). 



According to (|15|) and (|13|), in the case of standard viscoelastic solid model, the rebound 

(16) 



strain can be represented as 

er{t) =no- Ki exp(-(t - tm)/Ts 



where Ki is a dimensionless constant, which is proportional to the factor wo/h and depends 
on the ratios p and Tg/tm- 



The experimental points depicted in Fig. |5jwere obtained by Brown et al. (2009) in the 
indentation tests with a relative large maximum indentation strain of 30 per cent. Based 

(17) 



on (16), we make use of the following formula 
K2 - i^iexp(-(t - tm)/rs) 



Brit) 



with three fitting parameters Ki, K2, and r^. The results of the least-squares best fit 
of the analytical expression (17) to the experimental data are as follows: Ki = 0.116, 
Ts = 15.87 s, K2 = 0.234 (normal cartilage) and Ki = 0.107, r, = 15.65 s, K2 = 0.202 
(artificially degraded cartilage). It should be noted that in the both cases, the values of 
parameter K2 are different from 0.3 (which is the experimentally imposed maximum strain 
uq). In other words, the cartilage samples do not show complete recovery. Further, the 
characteristic relaxation time is found to decrease with degradation. Such a tendency 
is in agreement with that obtained by Hayes and Mockros ( 1971 ) in the shear creep test 
for human articular cartilage. 



6 Conclusion 



In the present study, the rebound indentation test has been evaluated in the framework 
of viscoelastic layer model. For the contact force P^^\t) in the displacement controlled 
stage and for the indentation displacement w^'^\t) in the load-controlled stage, the closed- 
form analytical expressions ([T]) and (10), respectively, have been derived for an arbitrary 



viscoelastic solid model. 
In the case of standard viscoelastic solid model, a simple exponential behavior was es- 



tablished for the rebound displacement w^^^{t) (see Eq. (13)). The theoretically predicted 
dependence of w^'^\t) on the time variable t only through the exponent exp{—{t — tm)/Ts) 



was verified using the three-parameter fitting formula for the rebound strain (17). 



From the view point of assessing viability of articular cartilage, the main finding of the 
present study is that the rebound displacement w^'^\t) does not depend on the layer 
material constants Eoo and u as well as on the layer's thickness h. At the same time, the 
rebound deformation mainly depends on the characteristic relaxation time r^, which is 
attributed to the articular cartilage permeability. 
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